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We compute the leading behavior of the uniform magnetic susceptibility, x, of a Fermi liquid near 
an antiferromagnetic transition with dynamic exponent z = 2. Our calculation clarifies the role of 
triangular "anomaly" graphs in the theory and justifies the effective action used in previous work . 
We find that at the z = 2 critical point of a two dimensional material, H?7ig^ox(3, 0) = Xo — DT 
with xo and D nonuniversal constants. For reasonable band structures we find that in a weak 
coupling approximation D is small and positive. Our result suggests that the behavior observed 
in the quantum critical regime of underdoped high-Tc superconductors are difficult to explain in a 
z = 2 theory. 



The anomalous magnetic behavior of the underdoped high Tc materials has kindled interest in the properties of 
antiferromagnetic quantum critical points. The data are most clear in Y Ba2Cu4^0s,- In this compound the planar 
copper nuclear relaxation rate If'^TiT ex 1/T^ for 2QQK < T < 6QQK while the uniform static susceptibility 
x(T) = Xo + DT^ in the same temperature range. The T-dependence of the copper 1/TiT in high T^ materials has 
been interpreted'* as due to strong, temperature dependent antiferromagnetic correlations, as expected if the spin 
dynamics is controlled by proximity to a T = antiferromagnetic critical point. Two possibilities have been proposed 
for the critical point; these are distinguished by the value of the dynamic exponent, z. One class of models is based on 
the z = 2 critical point expected for an antiferromagnetic instability of a fermi liquid system^. The value z = 2 arises 
because in a Fermi liquid the dominant contribution to the dynamics of an antiferromagnetic spin fluctuation is decay 
into a particle-hole pair. The z = 2 theory is argued**'® to apply to the high T^ materials because they are electrical 
conductors and presumably have fermi surfaces and therefore particle-hole continua into which spin fluctuations may 
decay. An alternative point of view^'® is that the correct model involves undamped (or weakly damped) spin waves 
with dispersion relation w = cfc, with the spin- wave velocity c remaining constant near the critical point, implying a 
dynamical exponent z = 1. For this point of view to be correct it is necessary either that the particle-hole continuum 
does not exist or that it does not couple to the spin waves. 

One phenomenological argument for the applicability of the z = 1 theory to Y Ba2Cu4^0s, involves the temperature 
dependence of the uniform susceptibility, x- Indeed, the result x(^) ~ ^^ has been demonstrated theoretically for 
the quantum critical regime of the z = 1 model®, and the coefficient A was shown to be universal (after the spin 
wave velocity and the gr-factor are fixed). The rough (factor-of-two) agreement between the predicted value oi A and 
the experimental results for dx/dT^ has led several authors '® to propose that the experimental susceptibility should 
be viewed as the sum of two terms: a weakly temperature dependent background term coming from a fermi liquid 
like continuum and a term proportional to T coming from quantum critical z = 1 dynamics. This point of view has, 
however, been questioned on phenomenological grounds®. 

In this communication we give a calculation of the critical behavior of x near a z = 2 antiferromagnetic critical 
point. In addition we clarify two related theoretical issues: the role of triangular "anomaly" graphs in the theory 
of the z = 2 transition and the role of "gauge invariance" arguments relating the response of a system in a uniform 
magnetic field to a stiffness against changes in boundary condition in the time direction-*-". 

Our theoretical approach follows previous work-*-'^: we integrate out the conduction electron degrees of freedom 
at fixed values of the staggered magnetization (j) and externally applied field hext, obtaining an expression for the 
effective action S^ff, as a functional of ^ and he^t- We then evaluate the free energy implied by S^ff using standard 
diagrammatic and renormalization group means. The result obtained in previous work at h^xt = is 

s'" = \Y. x-\k, n)${Q + k,n). ${-Q - k, -0] 

k,n 

+ U J2 <5(fcl+fc2+fc3 + fc4).5(Oi+02+03 + 04) (1) 

X [$iQ + fci, Oi) • ^{-Q + k2, n2)MQ + fcs, O3) • ^i-Q + h, O4)] + ... 

where x""*" is the inverse of the staggered susceptibility and is given by 

x-\k,n)=^-f + k' + 6 (2) 

Here F is a relaxation time, 6 measures the distance to the antiferromagnetic critical point, wavevectors are measured 
with respect to the ordering wavevector Q, and both the (p'^ coupling U and higher couplings are irrelevant (or 
marginally irrelevant) under renormalization in d > 2 spatial dimensions. 

The diagrams contributing to Seff at hext = are shown in Figs, la and lb. Diagrams with an odd number 
of (j) vertices are forbidden by momentum conservation if the momenta of the (j) fields are restricted to be near the 
ordering wavevector Q. The diagram in Fig. la is the sum of all diagrams which are irreducible with respect to the 
fermion-fermion interaction; this differs from the inverse susceptibility by an additive constant. 

The (fi'^ coefficient U is given by Fig. lb and is an analytic function of its arguments if the ordering wavevector 
Q is not an "extremal vector" of the Fermi surface. By "extremal vector" we mean one which connects two Fermi 
surface points with parallel tangents. For a circular Fermi surface, any Q = 2pF would be an extremal vector. The 
case where Q is an extremal vector will be discussed else where ■'-■'-. 

We now derive Seff in a nearly uniform external applied field. Before presenting the explicit calculations we recall 
arguments based on rotational invariance, which provide a useful check on our results. The essential point is that 
spins precess in an applied field, so the equation of motion of a spin Si in a field h is 



dSi/dt = h(Ri, t) X Si + i[H, S] (3) 

where H is the Hamiltonian in zero field. If the applied field is spatially uniform (but with an arbitrary time 
dependence) the all spins precess at the same rate, so if H is rotationally invariant the explicit dependence on field 
in the equations of motion may be removed by the rotation S"{t) = R^ (t)Si{t) with 



RTit) 



Texp dte'''"'h''{Ri,t') 
Jo 



(4) 



Therefore the time evolution operator in the presence of the field h, Uh(t), is given by (note R/j(0) = 1) 

Uh(t) = Kh(t)Uh=o(t) (5) 

where Uh=o{t) = e*^* . 2, From Eq. (5) it follows that the dependence of the time evolution on a spatially uniform field 
is only via boundary conditions. This line of reasoning applies also to the imaginary time formalism used to calculate 
thermodynamic quantities. Specifically, using Eq. (5) and the usual relation-*-^ between U(t) and the partition function 
Z we see that 

z(h) = j2<M^hmuh=om\n> (6) 

n 

Here n labels states in Hilbert space. If /i(t) = /i^e*""" with v = 2imiT ^ and integer m then from Eq. (5) we see 
Rj^(il3) = 1 and so we conclude from Eq. (6) that in this case d^Z/dh^ = 0. This will serve as a consistency check on 
our subsequent results. 

Note that at j^ = (i.e. for static h), Rh{i(3) ^ 1. However, in Appendix A we show by explicit calculation that for 
the model defined by Eq. (1), d^ Z/dhg = 0. This differs from the result found for the nonlinear sigma model in two 
dimensions, where d^ Z / OHq was shown to coincide with the physical susceptibility /ir77,g_>o[/i?77,tj_>ox(9) i^)]- 

Equation 6 can also be used to derive the dependence on magnetic field of the effective action S^ff describing the 
physics of low energy degrees of freedom if S^ff is local in time. In this case U can be represented as a Trotter 
product-*-^ and the quantity R}i{t)R^ (t) can be inserted at each time step yielding: 

N N 

Unm = n Riirj)e''^^R-\iir, - Ar)) = JJ e^^"^^^) (7) 

i=i i=i 

For example, if the low energy theory is the nonlinear sigma model with action (drn)^ + (\/n)^ , then it can be verified 
from Eq. (4) and Eq. (7) that this procedure leads to the replacement (S^n) —^ drn + h x n in the action, as it should. 

We have not been able to extend the arguments summarized in the previous paragraph to the z = 2 critical theory 
defined in Eq. 1. The essential difficulty arises from the fact that we have integrated out gapless fermionic degrees 
of freedom to obtain Seff- One consequence is that, as we shall see from the explicit calculations presented below, 
the coefficient coupling the field hextil,^) to the fields ^ is a nonanalytic function of q and as g, ^ and in 
particular depends strongly on the ratio 0/g. Such dependence is common in Fermi liquid theory and means that 
formulae such as Eq. (6), which are derived assuming spatially uniform fields with an arbitrary time dependence (i.e. 
0/g —f c»), may not yield the physical susceptibility, which corresponds to the opposite limit q —f with 0/g = 0. 

We turn now to an explicit calculation of S'^-^-^ in the presence of an external field. Because the four spin-wave 
interaction is irrelevant, only two diagrams must be considered: these are generated by inserting one or two hext 
vertices into Fig. la. The resulting diagrams are shown in Fig. 2. The analytic expression corresponding to these 
diagrams may be written: 



S'"{Kxt{q,in)) - S'"{Q) = iCuA(l^i^yabcKxM,i^)4>\Q + k,iv)4>'{-Q-k-q-iv-i^) 

D^^^iq, in) + iS"' - ^)i?W(g, iO)] h"Uq, in)h"U-q - iO) (8) 

b''(Q + k,iiy)<p''(-Q-k,-iiy) 



Here we have made the vector structure explicit and we have assumed that the fermion-fermion interaction producing 
the vertex corrections shown as the shaded areas in Fig. 2 is spin-independent. The coefficients C and D^^'^' are 
obtained by evaluating the diagrams shown in Fig. (2a) and Figs (2b, c) respectively. In the terms coming from Figs. 
(2b, c) we have retained only the contribution in which the two h^xt fields and two (j) fields have the same momentum 
and the same direction, because this is the only term we shall need in our subsequent calculations. 



2, From Eq. (8) we may evaluate x- Three diagrams arise; these are shown in Fig. 3. The resulting analytic 
expression is 

Xu(q, 0) = 2TJ2CI(1, n)x'(k, z.) - T^[3i?fj(g, 0) - D^,:^(q, n)]x(k, z.) (9) 

k,i/ k,i/ 

The sign of the term involving C^ ^(q, 0) has two contributions: i^ from the definition of the coupling constant and 
— 1 from the cross product. In order to evaluate Eq. (9) we require expressions for the coefficients C and D. We will 
use general arguments based on the cancellation implied by Eq. (6) for g = 0, 7^ 0. We then use Fermi-liquid 
arguments to obtain the relation of this limit to the limit = 0, g 7^ 0. We illustrate the general arguments with 
explicit calculations done in the weak coupling limit where all vertex corrections may be neglected. 
The general expression for C is 

Ck.iq,in) = 2Tj2^iP+<l,P)Gip,i€) [G{p+q,i€ + in)Ap,,{k,i^;q,n) (10) 

+ GiP + q,ie - ifl}Ap^eik, -v] q, -0)] 
The vertex correction F and block A are defined in Fig. 2a. In the weak coupling limit, T(p+ q,p) = 1 and 

Ap,,(fc,jy;g,0)= — -— -— (11) 

te + til -^iv - ep+Q+k+q 



2, From Eq. (2) we find 



Ck,„{^, in) = 2[x~\k, iv + iO) - x~'^{k, iv)]/iQ (12) 



CU0,iO) = ^^^?^ + i^Q" (13) 



with C''"(k, v) an analytic function of v and k^. 

The usual Ward identity arguments-'-'* imply that Eq. (12) is in fact correct as ^ for any model of fermions with 
a short range spin-independent interaction. The point is that diagrams for C are generated by inserting a magnetic 
field vertex in diagrams for x~^- At g = the magnetic field only shifts the energy of an electron line. The shift 
in energy may be absorbed into a shift in frequency and thus into a shift in external frequency, so a derivative with 
respect to magnetic field is equal to a derivative with respect to frequency. 

A similar argument shows that D^^> and D^''> can be expressed in terms of the second frequency derivative oi x~^- 
Consider for example a diagram giving the leading singular contribution to x""*" shown in Fig 4a. To obtain D^^> we 
assume that the external frequency is carried by the upper electron Green function and we differentiate twice with 
respect to the external frequency finding D^^> = 2— g^j— . To obtain D^''> we assume that external frequency is equally 
shared between the upper and the lower lines. These considerations imply that D^''> = —2 „^; . The numerical 
factors in these equalities arise as follows: D^ is defined to be the sum of the four diagrams of the type shown in Fig. 
2b, while the second derivative of x~^ produces two such diagrams if the external frequency is carried on the upper 
line; D^''> is the sum of two diagrams of the form shown in Fig. 2c, while the second derivative oi x~^ produces four 
diagrams of the type shown in Fig. 2b (if both derivatives act on the same line) and two of the type shown in Fig. 2c 
(if one derivative acts on an upper line and one on a lower line). The diagrams of the latter type enter with negative 
sign (because the external frequency flows backwards on the lower line) and there is a factor of 1/4 multiplying all 
diagrams (because one half the external frequency flows on each line). Combining all factors we obtain: 

3i?f )(0, 0) - i?g(0, 0) = 8[x-Hk, iy + iO) - 2x-\k, iv) + x'^ik, iv - iO)]/(iO)2 (14) 

or, using Eq. (2), 

3^^1.(0, 0) - i.i:)(0, 0) = 8 l- + »|-^H+l-»l + ,,«»(,, .) (15) 

By substituting Eqs. (12) and (14) for C and D into Eq. (9) and using the discrete analogue of integration by parts, 
one sees immediately that x(0, 0) = 0, consistent with Eq. (6). 



Because it will be important in subsequent arguments, we now discuss the origin of the nonanalytic frequency 
dependence in more detail. Consider first the term \v\lT in x~^- The usual Landau threshold arguments-*-^ imply 
that the nonanalytic behavior of x~^ comes from diagrams such as that shown in Fig. (4a), which may be divided 
into two disconnected parts by cutting two electron lines. The usual results of Fermi liquid theory imply that the 
vertex Vq coupling an antiferromagnetic spin fluctuation to the fermions and defined in Fig. 4 is an analytic function 
of fc, 0. The nonanalytic term |j^|/r comes from the product G(p, e) G{p-\- Q -\-k,e-\- v) and in particular comes from 
the region of the p, e integrals where both p and p + Q are near the Fermi surface and e is small. If more lines were 
present to share the frequency v, the nonanalyticity would be washed out. The general Ward identity says that in the 
0-limit inserting a field vertex is equivalent to differentiating with respect to frequency. Therefore, the nonanalytic 
contributions to C and D come from putting hext vertices on the isolated electron lines in Fig. (4a), as shown in Figs. 
(4b-d). 

We now turn to the limit g :^ 0, iO = which corresponds to the physically measured bulk susceptibility. Instead 
of calculating C(q, 0) directly we combine our previous result for C(0, 0) with a calculation of 

AC = C(g,0)-C(0,O) (16) 

Because the entire nonanalytic contribution to C(0, 0) is given by diagrams of the form shown in Fig. (4b), it is 
convenient to write 

C(q,n) = C4t,(q,n) + Crest(q,n) (17) 

with C46 given by the sum of all diagrams of the form shown in Fig. (4b) and Crest given by the sum of all other 
diagrams. We now calculate separately the contributions AC46 and ACrest- We begin with AC46, which is 

AC46 = TJ2[G{p + q, e)G(p, e)T(p +q,p)- G(p, e + 0)G(p, e)T(p + q, p)] (18) 

• TlG(p+Q + k,e+iy) 

Because the factor in brackets is convergent we may perform the p-integral first. In this integral the factor G(p + 
q, e)G{p, e) appearing in the first term contributes only for momentap far from the Fermi surface; for such momenta the 
v, k dependence is smooth. Therefore the first term in brackets gives only an analytic contribution to AC46. The second 
term in brackets is precisely 6*46(0,0). Therefore the nonanalyticity of AC46 precisely cancels the nonanalyticity in 
6*46(0,0) so that C*46(g, 0) is analytic. C*rest(g,0) is also analytic, for the same reason that C*rest(0,O) is analytic. 
This completes the argument that C*(g, 0) is an analytic function of v. Because it must be odd under time reversal 
we have 

AC*"" = iCiz^ + ... (19) 

Here C*i is a non-universal constant set by short length scale physics and the ellipsis denotes terms of order vk^ and 

Very similar arguments apply to D. The leading singular nonanalytic contribution to -D(0, 0) is cancelled by terms 
in AZ), so D{q, 0) ~ constant. There is one subtlety: among the diagrams in Drest are those with one field insertion 
on an isolated G-line in x~^ and one in Tq. In the 0-limit, the insertion in Tq gives a contribution dTq/dv. But 
because Vq is an even analytic function of j^, this gives a factor v, and therefore a contribution of order \i'\ to D. This 
is however negligible relative to the leading term, which is a constant. We therefore write ZD^ \q, 0) — D^ \q, 0) = D. 
Substituting into Eq. (9) one finds 

lim.^oxiq, 0) = 2TY,[iCiyfx\k, v) -Tj^Dxik, ly) (20) 

k,!/ k,i/ 

Using Eq. (1) in Eq. (20) we find 



^(g^0,0) = -2C,^T^ "' .2 -^E |.U^2^, +- (21) 



r 



Eq. (21) may be evaluated in different limits. One particularly important limit is the "quantum critical regime" in 
which (5 ^ T so that the critical behavior is cut off only by thermal fluctuations. It has been claimed'*'® that the 



underdoped high Tc superconductors are in this regime. Evaluating Eq. (21) in two dimensions, assuming 6 ^ T 
yields, to logarithmic accuracy, 

x{q^0,0) = Xbg + C2S- ^THT/6] (22) 

Here the positive constant C2 is determined by high energies where even Eq. (2) is not valid. 

In the quantum critical regime of the d = 2, z = 2 theory, it has been shown^ that 6 ~ T/ln[T] so that the 
dependence of x on T is, for all practical purposes, linear. 

In the d = 3 quantum critical regime the 6 dependence is not important because the k integral is convergent even at 
j^ = 0. Evaluating Eq. (21) by scaling one finds xil ^ 0, 0) = Xo + cx.DT^'^ + ... where a is a constant not determined 
by the scaling argument. 

We now discuss the sign and magnitude of D(q,0). We have calculated D by evaluating Figs. (2b) and (2c) for 
noninteracting fermions with spectrum €p and spin fluctuation vertex Tq. We obtain, after performing the frequency 
integration, 

D('')(Q 0) = 8 / '^'P f^^P^ - f(^P+Q) + (^P+Q - ^)f'(^p) + l(gg+9 - ep)V"(ep) (33) 

J (27r)2 {€p -€p+Qf 

i)(^){Q, 0) = -4 / '^'P 2(/(e,) - f{ep+Q)) + (e,+Q - e,)(f (e,) + f (e.+g)) ^^4) 

Here D^^''^) = D^^'''> /Tq. We have evaluated these integrals for tight binding bandstructures with dispersion 

€p = —2t(cospx +cospy) — it' cos px cos Py (25) 

and various ratios oit'/t. Some representative results are shown in Fig 5. We find that D = 3D^^>(Q, 0) — D^^'{Q, 0) 
is positive for |Q| < 2pf{0), negative for |Q| > 2pf{0). For \Q\ differing by more than Ri 20% from the critical value 
2pf{0) we find Z)/(47r) is of order 0.03rQ/t^. In a Fermi liquid one expects Tq to be set by an electron-electron 
interaction energy which is presumably ~ t, so Z)/(47r) ^ 0.03/t. 

The sign of D can be understood from the following qualitative arguments. First, from Eqs. (23,24) we see that if 
Q > 2pj?, Z) < because the terms involving derivatives of the Fermi function vanish and the sign of the denominator 
is opposite to that of the numerator. Of course, Q > 2pj? corresponds to undamped spin waves, which contradicts our 
starting assumption. Suppose now Q < 2pF. Because a Fermi liquid has no scale other than 2pF, the box diagrams 
must be smooth functions of the external momenta if these momenta are less than 2pF. Therefore it is reasonable to 
suppose that the sign of the diagram in 2b (when two external lines are at g Ri 0) is the same as lb (when all four 
external lines are at Q). But the sign of the diagram in lb is positive, corresponding to repulsive spin-fluctuation 
interactions and the stability of the theory defined by Eq. (1). Of course this argument is merely qualitative. It may 
be possible to devise band structures for which D < for Q < 2pF. 

To summarize, we have shown that in a two-dimensional Fermi liquid which undergoes a zero temperature z = 2 
antiferromagnetic transition the uniform susceptibility in the quantum critical regime is constant plus a term which 
is approximately linear in T. The coefficient and the sign of the leading T-dependence is nonuniversal. This leaves 
open the possibility that in real materials the temperature dependence of xil ~^ 0, 0) is determined by the terms 
which are non-leading in the z = 2 theory; for example the terms proportional to C2 in Eq. (22) which are small by 
a factor of 6/T. A discussion of these terms requires physics at high energy scales which is beyond the scope of the 
z = 2 critical theory studied here. 

However, for band structures similar to those proposed for the high Tc superconductors and for wavevectors Q not 
too close to 2pF, the sign is negative and the magnitude is small. The behavior we find differs in several respects 
for the behavior oi x(<l —^ 0,0) near the z = 1 transition discussed previously®. In the z = 1 case the transition is 
between a magnetically ordered phase with linearly dispersing undamped spin waves and a singlet phase with a gap 
to all spin excitations. The background term Xbg vanishes at T = at the critical point and in the singlet phase, and 
the leading term in the quantum critical regime is xil ^ 0, 0) ~ AT with A a universal and positive number. Note 
that if the ordering wavevector were outside the particle-hole continuum (Q > 2pj?) the spin fluctuation contribution 
to xil ~^ 0) 0) would be given by the universal z = 1 results but there would in addition be a fermion contribution 
described by Eq. (22) but with D < 0. The resulting xil ~^ 0) 0) might be in qualitative agreement with high Tc 
data. 

To conclude, we comment on the relationship of our result to previous work of Sachdev and to data. Sachdev has 
combined the rotation arguments discussed above Eq. (6) with finite size scaling to argue that near a quantum critical 

6 



point the free energy in an SU(2) invariant system is a universal function fu(hext/T) and that the "Wilson ratio" 
(susceptibility divided by specific heat coefficient) takes a universal value-*-". Our results do not contradict those of 
Sachdev; rather, they imply that in the physical limit 0/g —f 0, fu(hext/T) = 0, a possibility allo-wed by the general 
arguments. The mathematical reason for the vanishing of /„ is the nonanalyticity in q, of the coefficients coupling 
hextil,^) to the ordering fields (j). In the limit 0/g —f c» -where the rotation arguments apply the couplings have 
the expected form. Ho-wever, in the limit 0/g —f -which is relevant to measured bulk susceptibilities, the leading 
(universal) terms in the couplings vanish. The physical reason underlying this mathematics -was stated previously^: 
the universal results are ultimately derived from the fact that a uniform magnetic field generates rotations of, e.g. the 
staggered magnetization. Ho-wever, in z = 2 problems the staggered magnetization decays into a particle-hole pair 
before it has time to precess significantly, and so arguments based on precession do not give the leading behavior of 
x(q —^ 0, 0). A similar vanishing of the universal function is found in a class of impurity models-*-®. 

Our results also justify the action Eq. (1) used in previous -work-*'^ to analyze z = 2 transitions. This action 
neglected e.g. triangle graphs such as Fig. (lb) in -which one vertex is the nearly uniform magnetization and f-wo are 
the staggered magnetization. Our analysis implies that if these terms are included in S'^-^-^ (eq. (1)) and then the 
uniform magnetization fluctuations are integrated out, the resulting contributions are nonleading. 

Concerning experiments, we note that the observed linear T-dependence of the uniform susceptibility in high-Tc 
superconductors has been argued to be evidence of z = 1 critical behavior. We see that -would be rather difficult to 
explain in a z = 2 model, because for reasonable band structures the coefficient of the leading T-dependent term is 
negative and of small magnitude if the ordering -wavevector is inside the particle-hole continuum; ho-wever, there are 
non-leading (by a factor 6T/T) terms of positive sign -which could explain the data. 

Note added: After this research -was completed we received a preprint from S. Sachdev, A. Chubukov and A. Sokol 
reporting the study of the crossover from z = 1 to z = 2 behavior-*^. They obtain some-what similar results for 
x(q —^ 0, 0) but they do not consider logarithms and so do not distinguish bet-ween leading and non-leading terms. 
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APPENDIX: APPENDIX 

In this Appendix we evaluate Eq. (6) at T > for a static uniform magnetic field hext = hz. Equation (6) may be 
-written as a coherent state path integral, 

Z{h)= V<P*V<Pexp-Ty2 I d.'^ I dT'K{T-T')<Pk{T)-<p*_,,{T') 
J j^ Jo Jo 

(Al) 



Jo 



The fields (j) obey the rotated boundary conditions 

np) = R''\p)<i>\^) (A2) 

The kernel K{t) is given by 

ir(r) = T^e-|H/r (A3) 

v<r 

Because we have chosen the field to be in the z direction, the functional integral over ^^, (j)*' is /i-independent and so 
may be ignored. To perform the integral over the other components we define ne-w fields -f/)^,-!/)^ by 

^l (t) = cosh{hT)<t>l (t) - i sinh(/iT)^^ (t) (A4) 

^lt{T) = isinh(/iT)^^(T) + cosh{hT)<t>l{T) (A5) 

7 



We similarly define ip*" and ■f/'*^ using the complex conjugate of Eqs. (A4,A5). The transformation (A4,A5) has unit 
determinant, so 'D(j)'D(j)* = X'-f/'X'-f/'* . 2, From Eqs. (A2, A4,A5) and Eq. (5) we see ipiP) = V'(O)) so the ip fields may be 
expanded in Matsubara frequencies in the usual way, ipiT) = T'^^ e*""". 

We now substitute Eq. (A4,A5) into Eq. (Al), expand the fields ip in Matsubara frequencies, use Eq. (A3) and 
perform the t integrals. The only /i-dependence is in the first (double integral) term of the argument of the exponential 
in Eq. (Al). Denoting this term by ^j. Sk we find 

S„{h) = T^ J2 ^(2-2cosh/3/i)[(V'"(z.i)V"(z.2) + V'^(j^i)V'*^(j^2))/(j^i,J^2,J^3;M 



with 



(V'"(^l)V'*^(^2) + V'"(^l)V'^(^2)ff(^l, ^2, ^3; h)] 

2[(J^1 + i^3){i^2 + 1^3) - h^] 



g{vi,V2,V3]h) = - 



h[vi + 1^2 + 21^3] 

{{l^l + l^3f + h^Kii^2 + l^3f + h^) 



(A6) 

(A7) 
(A8) 



The first term (involving F) can contribute to Z to order h^ only if J^i = J^2- It is convenient to distinguish the two 
cases —i>3 = 1/1 = 1/2 and —1/3 zjz v\ = V2. One finds: 



si{K) = f 1 - ^) E ^(^.^(-)^r(-) + ^f (-)C(-)) 



-'^^EE 



12 y ^^ r 



V v^i^-v 



r {v^v^)'^ 



{i,l{v)i>T{v)^rMrA^)) 



(A9) 



The second term (involving g) plainly vanishes if v\ = V2 = —1^3 and is O(h^) if vi ^ 1^2 ^ —V3- We may therefore 
write 
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;,From Eqs. (A1,A9,A10) we obtain 
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(All) 



By use of the identity 
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along with the results v\ = 2imT and 
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11 = 1 



(A12) 



(A13) 



one may show that the last term precisely cancels the first two, so Z(h) — Z(h = 0) = 0. 
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FIG. 1. Diagrams contributing to the effective action for spin fluctuations, Eq. 1. The wavy line denotes a spin fluctuation 
and the shaded regions denote fermion two and four point functions which are irreducible with respect to the fermion-fermion 
interaction. 
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Fig 2 



FIG. 2. Diagrams describing the coupling between a nearly uniform external magnetic field (heavy dot) and fluctuations 
of the staggered magnetization (wavy line). The solid line denotes the electron propagator G and the shaded regions denote 
dressing of the diagram by the fermion-fermion interaction. Each diagram is understood to be the sum of two diagrams with 
opposite direction of electron lines. In Fig. 2a we have picked out the two fermion-magnetic field vertex T{p + q,p) (lightly 
shaded triangle) and the two fermion-two spin fluctuation vertex (heavily shaded trapezoid), because these quantities are 
needed for the discussion in the text. We have not similarly decomposed Figs. 2b and 2c because detailed expressions for these 
diagrams are not needed. 




Fig 3 



FIG. 3. Diagrams used in computation of contribution of spin fluctuations (wavy line) to uniform susceptibility. The shaded 
triangle in Fig. 3a is the vertex shown in Fig. 2a; the squares in 3b, c are the vertices shown in Figs. 2b, c. 



10 




n n 



Fig 4 



FIG. 4. Diagrams contributing to nonanalytic terms in x,C, and D. The heavy dot denotes the nearly uniform external 
field, the wavy line denotes the spin fluctuations, and the solid lines the dressed electron propagators. The shaded triangle 
is the vertex r(p + g,p) coupling the dressed electron propagators to the external magnetic field and the solid triangle is the 
vertex Tq coupling the dressed electron propagators to the antiferromagnetic spin fluctuations. 
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Fig 5 



FIG. 5. Plot of coefficient D' 
Eq. (25), t = 1, t /t = —0.1 and two electron densities: 



D{Tx + qy,0)/{'iTTl) 



calculated for noninteracting electrons with the dispersion given in 
= 0.5 (curve (a)) and n = 0.7 (curve (b)). 
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